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Summary. We review the physical and cosmological consequences of two possible 
electromagnetic couplings to the dark sector: (i) a neutral lightest dark-matter 
particle (LDP) with nonzero electric and/or magnetic dipole moments and (ii) a 
charged next-to-lightest dark-matter particle (NLDP) which decays to a neutral 
LDP. For scenario (i) we find that a relatively light particle with mass between 
a few MeV and a few GeV and an electric or magnetic dipole as large as ~ 3 x 
10~^^e cm (roughly 1.6 x 10~^ /xs) satisfies experimental and observational bounds. 
In scenario (ii), we show that charged-particles decaying in the early Universe result 
in a suppression of the small-scale matter power spectrum on scales that enter the 
horizon prior to decay. This leads to either a cutoff in the matter power spectrum, 
or if the charged fraction is less than unity, an effect in the power spectrum that 
might resemble a running (scale-dependent) spectral index in small-scale data. 

1 Motivation 

The origin of the missing 'dark' matter in galaxies and clusters of galaxies has 
been an outstanding problem for over 70 years, since Zwicky's measurement of 
the masses of extragalactic systems PP. Recent cosmological observations not 
only tell us how much dark matter exists but also that it must be nonbaryonic 
12] — it is not one of the familiar elementary particles contained within the 
standard model of particle physics. Dark matter is a known unknown. We do 
not know what the underlying theory of dark matter is, what the detailed 
particle properties of it are, nor the particle spectrum of the dark sector. 

Promising candidates for the lightest dark-matter particle (LDP) — those 
that appear in minimal extensions of the standard model and are expected to 
have the required cosmological relic abundance — are a weakly-interacting 
massive particle (WIMP), such as the neutralino, the lightest mass eigen- 
state from the superposition of the supersymmetric partners of the U{1) and 
SU{2) neutral gauge bosons and of the neutral Higgs bosons P), or the 
axion 0]. There is a significant theoretical literature on the properties and 
phenomenology of these particles, and there are ongoing experimental efforts 
to detect these particles. 

There has also been a substantial phenomenological effort toward plac- 
ing model-independent limits on the possible interactions of the LDP. For 
instance, significant constraints have been made to dark-matter models with 
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strong interactions 5^ and self-interactions [S], and various models with un- 
stable particles have been investigated [7] . Electromagnetic interactions have 
also been considered, and models with stable charged dark matter have been 
ruled out ^ while there are strong constraints on millicharged dark-matter 
models 

Dark matter is so called because the coupling of it to photons is assumed 
to be nonexistent or very weak. Here we ask the question, "How dark is 
'dark' ?" , and review several recent investigations that consider the physical 
and cosmological constraints to and effects of the electromagnetic interac- 
tions of the LDP and the next-to-lightest dark-matter particle (NLDP). In 
particular, in Section |2we discuss the consequences of a neutral LDP with 
nonzero electric and/or magnetic dipole moments TU", and in Section 13 we 
discuss the cosmological effects of a charged NLDP which decays to a neutral 
LDP [TTIIT^ . 



2 Dark-Matter: Electric and Magnetic Dipole Moments 
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Fig. 1. The constraints to \m^, {D, M)] that come from present-day searches and 
experiments. Viable candidates must lie in the shaded region. The short-dashed 
"relic abundance" curve shows where the dark matter would have a cosmological 
relic abundance f2^h^ — 0.135, assuming standard freezeout of annihilations via 
the dipole coupling to 7 and no X'X asymmetry. 
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In this section we consider the possibiUty that the dark matter possesses 
an electric or magnetic dipole moment. The result of Ref. ^Uli illustrated in 
Fig. m is that a Dirac particle with an electric or magnetic dipole moment 
of order ~ 10~^^e cm with a mass between an MeV and a few GeV can 
provide the dark matter while satisfying all experimental and observational 
constraints.^ 

The effective Lagrangian for coupling of a Dirac fermion x with a magnetic 
dipole moment M. and an electric dipole moment T) to the electromagnetic 
field is 

C^x = -\x<^^^'^iM+l,V)xF^r (1) 

Below we summarize various physical and cosmological limits to the form of 
interaction shown in Eq. (^). For further details of these limits see Ref. |1(J) . 

2.1 Dark Matter: Annihilation and Relic Abundance 

We assume x particles exist in thermal equilibrium in the early Universe and 
their dipole interactions freeze out when T drops below my^. Their cosmo- 
logical rehc abundance is fi^K^ ~ 3.8 x 10'' {m^/mp) In (^A/^J\n /A where 

A = 0.038^/g7'7ip;m^(CTannt') (see, e.g., Eq. (5.47) in Ref. 53)- Here, is 
the effective number of relativistic degrees of freedom at the freezeout tem- 
perature Tf ~ m^l A. x~X pairs annihilate to either photons or charged pairs 
through the diagrams shown in Fig. [21 and cTg^anV — a^^^2-fV + cr^^^ff = 
{V"^ +M'^)m\/2n -f N^eaiV^ +M'^), where iVoft is the number of /-/ pairs 
with ruf < m^. If H^h'^ = 0.135 then {V^ + M^)^/^ ~ 1.0 x lO'^^e cm for 
iriy, ^ 1 GeV, as shown in Fig. ^ The present-day mass density of x particles 
might differ from these estimates if other interactions are significant or there 
is a x~X asymmetry. 

2.2 Direct Detection 

In the nonrelativistic limit, the differential cross section for the process shown 
in Fig. is da/d[2 = Z^e^ {V'^ ^ M^) I{%'k'^v'^(\ ~ cos6')], where v is the 
relative velocity. Roughly speaking, a ^ [Zei^iT)^ -I- A^^)/27ru^ ~ 6.4 x 
10~^^ Z'^{'DIy + Ml'^) cm^, using v ^ 10^^ c.^ Current null searches in ger- 
manium detectors 14^ thus require (I?^? +-^17)^^^ ~ 10^'' at ^ 10 GeV 
— improving upon previous limits jl5j . 



^ We quote numbers for both the electric and magnetic dipole moments in units 
of e cm, where e is the electron charge. For reference, the Bohr magneton = 
eh/2me = 1.93 x 10~^^ e cm in these units. 

' [©17, Mit] = [V, >!]/(10-i^ e cm) 
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Fig. 2. Feynman diagrams for (a) annihilation of a x~X P^-"' to two photons and 
(b) x~X annihilation to charged /-/ pairs. 
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Fig. 3. Feynman diagrams for (a) the scattering of x particles by charged particles 
and (b) 7-x scattering. 



However, for large enough dipole moments, x particles will lose energy 
in the rock/shielding above the detector and evade detection in under- 
ground experiments. Given a shielding thickness L (in meters water equiv- 
alent), we obtain the lower bound V"^ + M.^ > [^m^w^ — ■j^Tj^pT^f^-E'th] x 

[^Lj:j.Zf li^^^^il + ,f^)]-\ where /iK,m] = m^m(m^ + m)-^ 



2m^ 

is the reduced mass, is the mass of detector nuclei, Eth is the threshold 
nuclear-recoil energy, the index i sums over the composition of the shielding 
material, and fi is the fractional composition by weight. The most restrictive 
constraints for large dipoles actually come from shallow experiments with null 
results such as the Stanford Underground Facility run of the Cryogenic Dark 
Matter Search JHj, and the Cryogenic Rare Event Search with Supercon- 
ducting Thermometers |17j . Airborne experiments, in particular the balloon 
experiment of Ref. ^H] and the rocket experiment of Ref. ^Hl provide im- 
portant complimentary constraints. The bounds due to all these experiments 
are shown in Fig. ^ 
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2.3 Constraints from Precision Measurements 

In Ref. Til the effect of the dipole interaction on the anomalous mag- 
netic moment of the muon, standard model EDMs, and corrections to Z- 
pole observables were considered. The strongest constraint was found to 
arise from the contribution of x particles to the running of a. Such run- 
ning affects the relationship between Gp, rnw, and the value of a at 
zero momentum: m^, — {TTa)/{\^GF) [{^ — rn'^Y/m'^){l — Ar)\ . In the 
standard model Ar^'^' = 0.0355 ± 0.0019 ± 0.0002 while experimentally 
^^Exp ^ 0.0326 ± 0.0023 yielding zir^'^"' < 0.003 with 95% confidence. 
The dipole interaction contributes to Ar via the diagram in Fig. ^ and so 
(V"^ + M^y^^ £ 3 X 10^^^ e cm is required, as shown in Fig.^ 




Fig. 4. One-loop correction to the photon self-energy induced by dipole moments 
M,T> oi the dark-matter particle. 



2.4 Direct Production 

Missing-energy searches for light (m^ ^ 1 GcV) dark matter in rare i?+ 
decays was suggested in Ref. |^ where Br{B^ iir++invisible) ^ 10"'' was 
derived. This limit requires P ^ 3.8 x 10~^^ e cm for < {mg+—mx+)/2 — 
2.38 GeV. Similarly, rare K'^ decays lead to the limit V ^ 1.5 x 10^^*^ e cm for 
< (mx+ — '7i^+)/2 = 0.18 GeV. These constraints are not yet competitive 
with other constraints shown in Fig. ^ 

In order to limit dipole couplings using collider experiments an expression 
for the rate // Xxx^ where X is some set of visible final-state particles, is 
necessary. Naive application of the effective Lagrangian in Eq. is invalid 
because perturbation theory breaks down when the energy scale of a process 
satisfies £ ^ l/V. Missing-energy searches from L3 {£ « 200 GeV) and 
CDF {£ = 1.8 TeV) can not be directly apphed to effective dipole moments 
T> > 10^^^ e cm and "D > 10"^"^ e cm, respectively, unless a high-energy 
theory is specified. 

2.5 Constraints from Large-Scale Structure and the CMB 

A dipole moment induces a coupling of the dark matter to the primordial 
plasma by scattering via the diagrams shown in Fig. |21 Dark matter couples 
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to the plasma at early times, and subsequently decouples. When coupled 
to the plasma, short-wavelength modes of the dark-matter density field will 
grow less quickly relative to the standard case. The long-wavelength modes 
that enter the horizon after dark matter decoupling remain unaffected. The 
full calculation in cosmological perturbation theory is provided in Ref . jlO| . 




multipole / 
(b) 

Fig. 5. The (a) matter (b) CMB power spectra including baryon-x drag for = 
1 GeV. The solid curve is for (D^ -|- = 1.4 x 10"^'^ e cm, short-dashed 

for (1)2 + >l2^i/2 = ;^ X 10-^*^ e cm, and long-dashed for {V^ + M^Y'"^ = 5 x 
10"^'' e cm. The data points are from (a) SDSS and (b) WMAP [21]. 
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In Fig.|Slwe show the hnear matter power spectrum and in angular power 
spectrum of the cosmic microwave background (CMB) for several values of 
the dipole moment and = 1 GeV. Using a Markov chain Monte Carlo 
algorithm (see, e.g., Ref. and data from SDSS ,22,, WMAP 23,, CBI 

VSA |2S], and Type la supernovae we find the bound shown in Fig. 
n Dipole moments as large asV ^ 10"^"^ e cm are thus cosmologically viable. 

2.6 Gamma Rays 

In the Galactic halo x~X pairs can annihilate to two photons through the 
diagrams shown in Fig. The non-observation of a gamma-ray line by 
EGRET leads to the limit shown in Fig. ^ A detailed search for a line flux 
with GLAST may find an observable signature for uiy^ ~ 0.1 — 1 GeV and 
Pit ^ 10. There may also be constraints to the annihilation rate from the 
low excess heat of Uranus ^7] . 

3 The Effect of a Quasistable Charged-NLDP 

In this section we discuss how the decay of a quasistable charged particle (f) 
(the NLDP) to a neutral dark-matter particle x (the LDP) suppresses the 
linear power spectrum on scales smaller than the horizon during the decay 
epoch. Prior to decay, the charged NLDPs couple to and oscillate with the 
primordial plasma. After decay, the plasma is coupled only gravitationally to 
the LDP. If all LDPs are produced through the late decay of charged NLDPs, 
then, as shown in Ref. JJ, the power spectrum is cutoff on small scales. 
For a lifetime r ~ 3.5 yr this would reduce the expected number of dwarf 
galaxies and may solve the small-scale structure problem of cold-dark-matter 
theory [2HI ■ If, instead, a fraction of LDPs are produced through the late 
decay of charged NLDPs, the linear power spectrum is suppressed only by a 
factor (1 — Z^)^. This suppression might be confused with a negative running 
of the spectral index as = dris/dhik in data that probes the power spectrum 
on small-scales. We describe these effects in further detail below. 

3.1 Charged Decay in the Primordial Plasma 

As — > X, the (f>/x comoving density drops/increases as p^a^ — m^n^^e"*/'^ 
and PxCi'^ = ^x^xoi^ ~ /^s"*^^) respectively. Here n^o = f^xPcrit I rn^ is the 
comoving density of dark matter, the fraction produced through 4> decays, 
a is the scale factor, and t is the cosmic time. The charged particles are 
tightly coupled to the baryons through Coulomb scattering. We can thus 
describe the combined 0-baryon fluid as a generalized baryon-like component 
p. 

In the synchronous gauge the perturbation evolution equations are iden- 
tical to the standard equations for the baryons (see, for example, Ref. j2Hp 



8 Kris Sigurdson 



with the subscript b replaced by (3. For the dark matter — ~ikV^ — ^/i + 
^m^Jih-5x)^ and = ^^^v^ + Xr.'fJ[Vp -V^), where 5^ = 5pjp^, 
is the bulk velocity, Am = m^/m^^ and an overdot is a derivative with respect 
to conformal time. The modifications to photon perturbation evolution are 
negligibly small. 

Due to Compton scattering the (3 component and the photons are tightly 
coupled as a /3-photon fiuid at early times which supports acoustic oscilla- 
tions. Since dark-matter perturbations are sourced by the perturbations of 
the (3 component, /c-modes that enter the horizon prior to decay, when p^j 
is large, will track the oscillations in the /3-photon fluid rather than growing 
due to gravity. After decay, when p^j Px is small, dark-matter modes that en- 
ter the horizon undergo the standard growing evolution. In Fig. Elwe follow 
the evolution of the dark-matter perturbations through the epoch of decay. 
We used a modified version of cmbf ast 30 . 
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Fig. 6. The comoving k = 30.0 Mpc~^ (left), k = 3.0 Mpc~^ (center), and 
k — 0.3 Mpc~^ (right) 5^ perturbation for the ylCDM model (long-dashed) and a 
f(P = 1 model with r = 3.5 yr (solid). Also shown is the Sfs perturbation (short- 
dashed) . 

In Fig. we plot the linear power spectrum for — 1 and r = 3.5 yr. 
Power is cutoff for k~-^ ^ 0.3 Mpc relative to a standard ylCDM power spec- 
trum, reducing the expected number of subgalactic halos and bringing pre- 
dictions in line with observation In Fig. |7|3 we show the linear power 
spectrum for several values of r and J^. The linear power spectrum is now 
suppressed by a factor of (1 — f^Y- On large scales the linear power spec- 
trum describes the statistics of density fluctuations, but on small scales the 
full nonlinear matter power spectrum is required. The effects of nonlinear 
evolution are accounted for in Ref. ^21 where it is shown that the effect of 
the charged-NLDP decay on the small-scale nonlinear power spectrum can be 



How Dark is 'Dark' ? Electromagnetic Interactions in the Dark Sector 9 

similar (but different in detail) to that of a model parametrized by a running 
of the spectral index ag = dng/dhik. 

We note that since the NLDP is charged, LDP production will be accom- 
panied by an electromagnetic cascade. The latter could in principle reprocess 
the light elements produced during big bang nucleosynthesis, or induce unrea- 
sonably large spectral distortions to the CMB. In fact, the models discussed 
are safely below current limits [111 IT^ . 

3.2 Particle-Theory Models of Quasistable Charged NLDPs 

In the — 1 case, in order to be a solution to the small-scale structure 
problem the charged-NLDP must have a comoving density equal to the dark- 
matter density today and r ^ yr. One way to make the NDLP long lived is 
to suppress the LDP coupling. If LDP is a stable super-weakly-interacting 
particle [22], such as a gravitino G or the first Kaluza-Klein graviton 
in the universal extra-dimensions scenario the NLDP can be charged 
and, since it decays gravitationally, have a super-weak decay rate. The re- 
quired mass-spliting and lifetime can be achieved for « ~ 100 TeV. 
These masses are above the TeV range discussed in most supersymmetric 
phenomenology, and uncomfortably close to violating the unitary bound for 
thermal production |84j — nonthermal production is likely required for such 
large masses. 

For < 1, remarkably, there are configurations in the minimal super- 
symmetric extension of the standard model (MSSM) with the properties re- 
quired here |12j . If the LDP is a neutralino quasi-degenerate in mass with the 
lightest stau, one can naturally obtain, at the same time, LDPs with a relic 
abundance fi^h'^ = 0.113 j2] and NLDP lifetimes and the densities needed 
in the proposed scenario. Such configurations arise even in minimal schemes, 
such as minimal supergravity (mSUGRA) |35| and the minimal anomaly- 
mediated supersymmetry-breaking (mAMSB) model A detailed study 
of the (t,/^) parameter space using current and future cosmological data 
may constrain otherwise viable regions of the MSSM parameter space. Fur- 
thermore, this scenario may be tested at future particle colliders (such as the 
Large Hadron Collider (LHC) ) or dark matter detection experiments (See 
Ref. ^21 for full details) . The decays of staus at rest might even be studied by 
trapping these particles in large water tanks placed outside of LHC detectors 
ISZ). Other methods of trapping staus have also been considered |SH1- 

4 Conclusion 

We have considered two distinct scenarios for electromagnetic interactions in 
the dark sector. In the first we considered the effects of a LDP with nonzero 
electric and/or magnetic dipole moments and found that a light particle 
with Ri MeV-GeV and a dipole as large as 3 x 10~^^e cm « 1.6 x 
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Fig. 7. (a): The linear power in the standard ylCDM model (dashed), and a = 1 
charged-NLDP model (solid) with r = 3.5 yr. (b): Shown is Zi^(fc) = k'-^P{k)/2n^ 
for a ria = 1 /ICDM model (dashed) and for several charged-NLDP models. 



10^^ fiB satisfies experimental and observational bounds. The second scenario 
examines how the decay of charged NLDPs in the early Universe results 
in, depending on f^, either a cutoff or suppression of the small-scale power 
spectrum. For a purely gravitationally interacting LDP configurations with 
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= 1 may exist, while configurations with f^<l can be found within the 
MSSM. 
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